Cosmic microwave background with Brans-Dicke gravity II: constraints with the 

WMAP and SDSS data 
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Using the covariant formalism developed in a companion paper (paper I), we derive observa- 
tional constraint on the Brans-Dicke model in a flat Fricdmann-Lemaitre- Robertson-Walker (FLRW) 
universe with cosmological constant and cold dark matter. The cosmic microwave background 
(CMB) observations we use include the Wilkinson Microwave Anisotropy Probe (WMAP) five year 
data, the Arcminute Cosmology Bolometer Array Receiver (ACBAR) 2007 data, the Cosmic Back- 
ground Imager (CBI) polarization data, and the Balloon Observations Of Millimetric Extragalactic 
Radiation and Geophysics (BOOMERanG) 2003 flight data. For the large scale structure (LSS) we 
use the matter power spectrum data measured with the luminous red galaxy (LRG) survey of the 
Sloan Digital Sky Survey (SDSS) Data Release 4 (DR4). We parameterize the Brans-Dicke param- 
eter u) with a new parameter = ln(l/aj -I- 1), and use the Markov-Chain Monte Carlo (MCMC) 
method to explore the parameter space. We flnd that using CMB data alone, one could place some 
constraint on positive C, or u, but negative C, or u could not be constrained effectively. However, with 
additional large scale structure data, one could break the degeneracy at < 0. The 2(j (95.5%) bound 
on C is -0.00837 < C < 0.01018 (corresponding to lj < -120.0 ov lu > 97.8). We also obtained con- 
straint on G/G, the rate of change of G at present, as -1.75 x 10"^^yr~^ < G/G < 1.05 x 10~^^yr~^, 
and 5G/G, the total variation of G since the epoch of recombination, as —0.083 < SG/G < 0.095 
at 2(7 confidence level. 



I. INTRODUCTION 

The Jordan-Fierz-Brans-Dicke theory [1 H, i, i, i] 
(hereafter the Brans-Dicke theory for simphcity) is the 
most natural alternative to the standard general relativ- 
ity theory and the simplest example of a scalar-tensor 
theory of gravity 0, i, i [S [HI- 

The gravitational con- 
stant becomes a function of space and time, and is pro- 
portional to the inverse of a scalar field. Its action in the 
usual (Jordan) frame is 



where T^™'' is the stress tensor for all matter except for 
the Brans-Dicke field, and the equation of motion for ip 
is 



^ rp(ni)fJ. 

2lu + 3 ^ ' 
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In order to match the result of Cavendish type experi- 
ments, the present day value of ip should be 



2uj- 



2cj + 3 



(5) 
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where is the Brans-Dicke field, w is a dimensionless 
parameter, and 5''™-' is the action for the ordinary matter 
fields 5^'"^ = / d'^x^P^C}-™'^ . For convenience, we also 
define a dimensionless field 



(2) 



where G is the Newtonian gravitational constant. The 
Einstein equations are then generalized to 



G^. = — rW-H4(V,,^V.(p-l.g^.VA^VV) 

ip (^^ z 
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The original motivation of the Brans-Dicke theory is 
the idea that the gravitational constant G ought to be 
related to the average value of a scalar field, which is 
determined by the mass density of universe, so that the 
Mach principle is satisfied 0,0]. Later, it is noted that 
the scalar-tensor gravity appears in the low-energy limit 
of supergravity theories from string theory and other 
higher-dimensional gravity theories p^ . The Brans- 
Dicke field may be associated with the dilaton-graviton 
sector of the string effective action [l^, [31 ■ The dimen- 
sionless parameters in string theory - including the value 
of the string coupling constant - can ultimately be traced 
back to the vacuum expectation values of scalar fields 

The unexpected discovery of the accelerating expan- 
sion of the Universe [H, [3 forced us to look for 
an explaination of the so called "dark energy" which 
may drive such acceleration. Scalar fields rolling down 
a proper potential may serve as a dynamical dark energy 
model JJi, 2Qi, 3lL 22, 23, 2± 2^. However, in these 
phenomenological models, the scalar fields are added by 
hand, the connection to fundamental physics is often 
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unclear. The Brans-Dicke field is a natural candi- 
date for the scalar field, this is the so called "extended 
quintessence" scenario jS, 113, S [H S lU, [3||. Al- 
ternatively, the Brans-Dicke theory could also serves as 
an effective model of the f{R) gravity, in which the 

Pvity is invoked to ex plai n the cosmic acceleration 
, m m, m, m m mm, \m\mm\mmm,\m- 

The Brans-Dicke theory is reduced to the Einstein the- 
ory in the limit of 



0, ^" 



0. 



(6) 



So in some sense it could never be excluded completely if 
Einstein's general relativity theory turns out to be the fi- 
nal words on the classical theory of gravitation. So far, no 
significant deviation from the Einstein theory has been 
discovered, and the most stringent limit on the Brans- 
Dicke theory comes from solar-system experiments which 
constrain the parameterized post-Newtonian (PPN) pa- 
rameter J = {1 + oj)/ {2 + uj). A recent significant result 
was reported in 2003 using the Doppler tracking data of 
the Cassini spacecraft while it was on its way to Saturn, 
with 7 - 1 = (2.1 ±2.3) x 10"^ at 2a confidence level [H, 
which corresponds to about |cj| > 40000. The limitation 
of such experiments is that they are "weak-field" experi- 
ments and probe only a very limited range of space and 
time. They could not reveal spatial or time variation of 
gravitational constant on larger scales. 

It has long been known that cosmological observations 
such as the cosmic microwave background (CMB) and 
large scale structure (LSS) could be used to test the 
Brans-Dicke theory [11, Si, [H HH IH, [H [H [H, M, 
[stI . [ssl [sol . While the constraints obtainable with such 
methods are generally weaker than the solar system tests, 
they probe a much larger range of space and time. In re- 
cent years, with the launch of the WMAP satellite, and 
the completion of the 2dF and SDSS redshift surveys, it 
is interesting to put such test into practice. 

In 2003, Nagata et al. used the WMAP first year 
data and test method to derive a constraint on the 
Brans-Dicke parameter. They obtained uj > 1000 at 2a 
confidence level[55|. However, in 2004, Acquaviva et al. 
obtained a new constraint using a Markov Chain Monte 
Carlo approach with CMB data from the WMAP first 
year data, the ACBAR, VSA and CBI data, and the 
galaxy power spectrum data from 2dF. They obtain a 
result of w > 120 at 2a confidence level[56j. These two 
limits differ by an order of magnitude. We are unable to 
reproduce the result of Ref. [55|, but we did reproduce 
successfully the result of Ref. [56] using the procedures 
described in their paper and the same data set as they 
used. 

Nevertheless, as will be discussed in the next section, 
there is room for improvement upon the method used 
in Ref. [s^. Moreover, new CMB and LSS data have 
since become available, it is therefore time to revisit this 
problem with a new approach and update the constraint 
with the latest observational data. 

We have developed a covariant and gauge invariant 



method for calculating the CMB anisotropy in Brans- 
Dicke theory, the formalism of our approach is presented 
in the companion paper [l| (paper I). In the present pa- 
per, we apply the method developed in paper I, and use 
the latest CMB data and large scale structure data to 
constrain the Brans-Dicke parameters. Here we consider 
only the case of the massless Brans-Dicke model with 
cold dark matter and cosmological constant. The more 
interesting case of the Brans-Dicke field with interacting 
potential would be investigated in future study. 



II. METHODS 

The formalism of calculating CMB angular power spec- 
tra and matter power spectrum in Brans-Dicke theory 
with the covariant and gauge invariant method are pre- 
sented in paper I. We also described in that paper the 
numerical implementation of the method in the CMB 
code CAMB 60]. The results of the modified CAMB 
code have been check with the results given by Chen 
and Kamionkowski (1999) in Ref.fsS], which was based 
on a modified version of CMBFAST in the synchronous 
gauge. The output of the two code show excellent agree- 
ment. Our new code has been implemented with some 
techniques to improve the architecture of program, and 
the code is much faster than the old one. We refer the 
reader to paper I for more details. 

We consider deriving the constraint on the Brans-Dicke 
model with the observational data using the Markov 
Chain Monte Carlo (MCMC) simulation. The CAMB 
code is used by the public COSMOMC code [iH as the 
driver for calculating the CMB angular power spectra 
and matter power spectrum. Here we use the modified 
CAMB code in the COSMOMC simulation. 

The data we used to constrain the Brans-Dicke 
model are the latest cosmic microwave background 
power spectrum data, which include the WMAP five- 
year '62], ACBAR 2007163], CBI polarization [H and 
BOOMERanC 2003[6l M, 113 data. We also use the 
galaxy clustering power spectrum data derived from the 
SDSS LRG survey DR4 68]. 

We do not use the Type la supernovae (SNe la) data 
when making the constraint in this paper, because the 
value of the gravitational constant varies during the 
expansion of the Universe. We know that the Chan- 
drasekhar mass Mch esc G~^l^ . The variation of the 
gravitational constant G means that the peak luminos- 
ity of SNe, which is approximatively proportional to the 
Chandrasekhar mass, will change, so the supernovae can 
not be assumed to be standard candles in this model. 

Besides the Brans-Dicke parameter, the cosmological 
parameters explored in our MCMC simulation are {r^b/i^, 
^rrJr?, e, T, Hs, log{10^° A^), Agz} ■ ^bh^ , ^mh"^ are the 
baryon and matter densities respectively. The 9 param- 
eter represents the ratio between the sound horizon and 
the angular diameter distance to the last scattering sur- 
face, it is used in lieu of the Hubble parameter h since 
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it is less correlated with other parameters, r is the op- 
tical depth to reionization, Ag is the amplitude of pri- 
mordial superhorizon power spectrum in the curvature 
perturbation on 0.05Mpc~^ scale, is the scalar spec- 
tral index, Asz characterizes the marginalization factor 
of Sunyaev-Zel'dovich effect. We only consider the Brans- 
Dicke model in a flat universe with the cosmological con- 
stant as dark energy. We assume flat priors for these 
parameters, and the allowed ranges of the parameters 
are wide enough such that further increasing the allowed 
ranges has no impact on the results. 

In any Bayesian approach to error estimate and param- 
eter constraint, the result will depend somewhat on the 
parametrization and prior. The original Brans-Dicke pa- 
rameter uj is inconvenient to use, because it is unbounded, 
and the Einstein limit appears at w — > cx). Even if one 
restrict the allowed range of uj to some finite interval, 
the large to region would be unduely favored, because at 
such region the difference in CMB and LSS produced by 
models of different to becomes indiscernibly small. 

Acquaviva et al. introduced a variable In^ — 
ln[l/(4w)] in Ref . . and set its prior to be uniform 
in the range In^ G [—9,3], corresponding to w e 
[0.01, 2025.77]. The choice of the lower limit of In^ is mo- 
tivated by the fact that for cu > 2000, visual inspections 
show that the CMB angular power spectra become in- 
sensitive to Lo. This parametrization is workable, but has 
some drawbacks: firstly, it does not include the negative 
values of w, and secondly, the lower limit of ln^mi„ = —9, 
while ostensibly a reasonable choice, is nonetheless put 
in by hand and quite arbitrary. In fact, the 2a limit 
would be sensitive to this artificial choice, because the 
likelihood is high and almost flat at Inf < —9, so if one 
varies the lower limit ln^„ii„, the overall normalization 
of the posterior probability distribution function would 
be directly affected. 

In this paper, we introduce a new parameter which is 
more convenient to use: 

C = ln(l + -). (7) 

This parameter has the nice property that C ^ asymp- 
totes the Einstein gravity, and it is easy to obtain the 
two-side (i.e. allows negative uj) likelihood distribution 
around C = 0. C~l/w when w is a large number (i.e. 
very close to Einstein gravity). We set the allowed range 
as S [—0.014,0.039], which brackets the Einstein grav- 
ity case, and corresponding to lu £ [— oo,— 71] U [25, oo]. 
There is no arbitrary limit on large \uj\ value, but only 
limit on small value. Outside this range, i.e. —71 < 
w < 25, our numerical code break down, because the 
background evolution deviates too much from the stan- 
dard model. However, as we are looking primarily for 
small departures from the Einstein gravity, this is not a 
big concern, and large departures would have been eas- 
ily detected by other means as well. When making plots 
of the likelihood, we do taken into account of the range 
of allowed parameters, so that the probability is prop- 
erly normalized. Unavoidably, this artificial restriction 



on parameter range has some effect on the final result, 
but as long as the final probability distribution is much 
smaller than the allowed range, it would not fundamen- 
tally change our conclusion. 




FIG. 1: The one dimensional marginalized likelihood distri- 
butions for the parameter (. "WMAP5" denotes WMAP 
five-year data. "ALL CMB" represents WMAP five year 
data plus some small scale data and polarization data, i.e. 
ACBAR 2OO7EII, CBI polarization 64] and BOOMERanG 
2003[63, l&sl . l67l| data. "LSS" means galaxy clustering power 
spectrum data from SDSS DR4 LRG data. 



III. RESULTS 
A. Constraint on Brans-Dicke Theory 

The one-dimensional marginalized likelihood distribu- 
tions for C is shown in Fig[TJ The three curves are ob- 
tained with the WMAP data alone (magenta dash-dot 
curve), with all CMB data, i.e. WMAP, ACBAR, CBI 
and Boomerang data (blue dashed curve), and with all 
CMB data as well as the LSS data from SDSS LRG sur- 
vey (red solid curve). Interestingly, using only the CMB 
data, we find that a negative C, is favored. Indeed, the 
two curves obtained with only the CMB data declines 
very slowly at C < 0, making it difficult to obtain a limit 
on negative C with them, so we could not easily quote 
a number for CMB-only constraint. However, with the 
additional constraint from the large scale structure data, 
the best fit value of C goes back to the neighborhood of 
zero, and the likelihood declines rapidly (almost Gaus- 
sian) at negative (. This shows that the large-scale struc- 
ture data play an very important role in constraining the 
Brans-Dicke gravity. 

To understand this result in more detail, we consider 
three models: 

(1) the original best fit minimal 6-parameter ACDM 
model with Einstein gravity obtained by the WMAP 
team using their 5 year CMB data combined with the 
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FIG. 2: CMB angular power spectra data with predictions of 
three best fit models, see text for details. 



WMAP ACDM b=1.9 
ALL CMB + LSS b= 1 .9 
WMAP5 b=1.9 
WMAP5 b=2.2 




Figl3] respectively. 

As shown in Figl2l due to parameter degeneracy, the 
differences between the three curves of CMB are almost 
indiscernible: for a slightly negative the Brans-Dicke 
model could produce CMB spectra which fits the data 
very well. However, as shown in Figl3l the matter power 
spectra are quite different. The Brans-Dicke model which 
best-fit the CMB data does not fit the galaxy power spec- 
tra very well. To be sure, if one also allow the bias param- 
eter as a free parameter, the fit could be somewhat im- 
proved, nevertheless, it still fails compared to the model 
obtained by fitting both the CMB and LSS data. Thus, 
we see that the galaxy power spectrum data could play 
an important role in distinguishing models, even though 
when used alone its constraining power is relatively weak. 

The 95% marginalized bound we derive in this paper 



0.00837 < C < 0.01018 



corresponding to 



< -120.0 or uj> 97.: 



(8) 



(9) 



We note that when comparing this result with that of 
Ref. [sHi, one must remember that we have adopted 
different parametrization and priors. In fact, we have 
used CMB data with higher precision (WMAP 5 year vs 
WMAP 3 year), and additionally we also used the LSS 
data (SDSS) which they did not use. Despite this im- 
provement in data quality, the limit we derived appears 
to be slightly weaker than theirs, this is due to the differ- 
ent parametrization and prior we used, particularly, we 
allowed negative uj which was not considered in Ref. [56| . 

To better understand the degeneracy and the 2-D like- 
lihood space distribution, we plot the 2-D contours of the 
marginalized likelihood distributions of C against Q\ in 
FigHl The Einstein gravity with JIa ~ 0.75 is still the 
best fit model for the all CMB+LSS data set. If C is 
greater, should also be greater, and vice versa. 



FIG. 3: The linear galaxy power spectra given by the three 
best-fit models compared with SDSS LRG DR4 datafsH]- We 
adopt original value of bias of b — 1.9 as given in Ref. [6^. 
For the best fit model using only WMAP 5 year data, we also 
plot the result adjusting b value to 2.2 to better fit the galaxy 
power spectrum for comparison 



distance measurement from SN and the Baryon Acous- 
tic Oscillations(BAO) in the distribution of galaxies [69|. 
which is marked as "WMAP ACDM" in the figure; (2) 
the best fit Brans-Dicke mode using only WMAP five- 
year CMB data, which is marked as "WMAP5" in the 
figure; (3) the best fit Brans-Dicke model using all CMB 
data as well as the SDSS LRG data, which is marked as 
"All CMB-hLSS" in the figure. 

The CMB angular power spectra and linear galaxy 
power spectra for these models are plotted in Fig[2] and 




ALL CMB+LSS 
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FIG. 4: The 2D contours of the marginalized likelihood dis- 
tribution of ^ against S7a- 
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B. Constraint on cosmological parameters 




FIG. 5: The one dimensional marginalized probability distri- 
butions for the other cosmological parameters in Brans-Dicke 
theory and in General Relativity. Data are all CMB data and 
the LSS data from SDSS LRG survey (i.e. ALL CMB+LSS). 
Red solid curves are results for Brans-Dicke gravity, black 
dotted curves represent the case for ACDM model in General 
relativity. 

In Figini we plot onc-dimcnsional marginalized likeli- 
hood distributions for other parameters in Brans-Dicke 
theory (red solid curves), for comparison, we also plot the 
same distributions in General Relativity case(black dot- 
ted curves) which fixes ^ = 0, using the same dataset — 
"ALL CMB-I-LSS" , i.e. all CMB data combined the LSS 
data from SDSS LRG survey. The parameters in the top 
two rows of panels are the primary cosmological param- 
eters used in the MCMC program, and the parameters 
in the bottom two rows of panels are the derived pa- 
rameters (not the parameters really run in the MCMC 
code). We see that the best fit value of the parame- 
ters are almost unchanged. Furthermore, for most of the 
primary parameters, the width of the likelihood distribu- 
tion is also unchanged. Only the distribution of the dark 
matter density parameter fic^^ is slightly broader. For 
the derived parameters, the best-fit values are also ba- 
sically unchanged. However, the likelihood distribution 
for most parameters are broader, showing the introduc- 
tion of the Brans-Dicke model allows larger uncertainty 
in these parameters. The notable exception is the reion- 
ization redshift z^e which is basically unaffected. 

The 2-D contours of of the marginalized likelihood dis- 
tributions of Q against other cosmological parameters are 
shown in Figl6l As can be seen in the upper two rows 
of panels, there are apparently not much correlations be- 
tween C and the other primary cosmological parameters 
used in MCMC program, such as ilfc/i^, i^mh^, t, Us 
and \og{10^° As) ■ However, from FiglHand the lower two 
rows of panels of FiglSl we see that ( is correlated with a 



and the derived parameters including the age of Universe, 
Hqi r^m, and CTg, though there is almost no correlation 
with the reionization redshift z^e- 

We summarize the 68% confidence limits on cosmolog- 
ical parameters in Table[Tl Note that our pivot wavenum- 
ber fco = 0.05 Mpc^^ of the primordial power spectrum 
is different from that of the WMAP group 5 year data 
release (fco = 0.002 Mpc~^), and the set of primary pa- 
rameters we used is also slightly different from the one 
used by the WMAP group, as they used instead of 
as a primary parameter. As we have mentioned, the 9 
parameter is less correlated with C, hence our choice in 
this case could help improve the efficiency of the MCMC 
method. The data used by the WMAP group [6§| are the 
WMAP five-year data. Type la supernovae data and the 
BAO data. We have not included the supernovae data, 
which we considered unreliable in the case of modified 
gravity. From the Table HI we find that our best-fit values 
of cosmological parameters are generally consistent with 
the WMAP group result at one a confidence level, how- 
ever, our constraints are a bit weaker than those given 
by the WMAP group, as we have added the Brans-Dicke 
parameter, and also used somewhat different data sets. 

C. Constraint on the variation of gravitational 
constant G 

An interesting question is what limit could one place 
on the variation of the gravitational constant G using the 
CMB and LSS observations. In the Brans-Dicke theory, 
G also underwent evolution from the time of recombina- 
tion to the present time, the variation in G is correlated 
with the value of C, so we can also derive a limit on the 
variation of the G. Of course, this evolution is not arbi- 
trary, but determined by the dynamical equation Eq. 
so when citing the bounds on variation of G obtained in 
this way, one has to remember its limitations. Neverthe- 
less, we note that in the Brans-Dicke theory, the impact 
on CMB and LSS comes mainly from the variation of G 
[l|, [5^ , so the result obtained this way could still serve 
as a good reference value. 

For making this constraint, we introduce two derived 
variables in the MCMC, namely the rate of change of 
the gravitational constant G/G at present and the inte- 
grated change of gravitational constant since the epoch 
of recombination SG/G: 

G/G = -^/^, SG/G = {Grec - Go)/Go. (10) 

The one dimensional marginalized likelihood distribu- 
tions of G/G and 6G/G are plotted in Fig. [7] and Fig. 1 
respectively. The "WMAP 5 year data" and the "all 
CMB data" both favor a slightly non-zero (positive) 
G/G. With the addition of the SDSS power spectrum 
data, however, the best-fit value is back to zero. From 
these figures, we could still see some effect of the prior, as 
the likelihoods are still non-zero or at best just approach- 
ing zero at the edge of the figures. Nevertheless, with the 
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FIG. 6: The 2D contours of the marginalized likelihood distribution of against the other cosmological parameters. The color 
is for the mean likelihood distribution. 



TABLE I; Summary of cosmological parameters and the corresponding 68% intervals. 



Class 


Parameter 


WMAP5 


ALL CMB 


ALL CMB+LSS 


WMAP group [69] 


Primary 




02190+" 


02200+" ™°''^ 


n 09990+0-00033 


0.02265 ± 0.00059 






U.lU'iU_Q 0043 


1064+°-™" 

U.lUO^__Q QQ3Q 


n 1 nc^ + 0.0042 
0.1066„o 0046 


0.1143 ± 0.0034 




e 

T 


1 0391+''-°°*^ 

ll.U8»_(, 


1 0425+°-°''^^ 

n os=;+°-''" 

U.UOO_Q 


1 0432+''-°°^'* 

l.U«iOZ_o 0030 

093+" '"'^ 

U.UMO_o 014 


0.084 ± 0.016 




Us 


047+0 035 




962+" "^^ 
u.aDz_o.oii 


960+°"" 
u.yuU o.oi3 






^04+0. 074 


J.UOU_Q Q24 


3 070+0 030 
J.UrU_o.047 


A, = (2.457t;; »^^) X 10-° 


Derived 




o.78o;s::™ 


n 7SQ + 0-076 


n 7CO + 0.029 


0.721 ± 0.015 
0.0462 ± 0.0015 
0.233 ± 0.013 




Age/Gyr 




13.821°;^^ 


13.63l°-f4 


13.73 ± 0.12 Gyr 




n„ 




n 91 0+0.093 


0.247t«o:Si;, 












0.789l«°o55 


0.817 ± 0.026 








10.2i;;» 


lo.gt";^ 


10.8 ± 1.4 




Ho 




64 4+"-^ 


72.3t4° 


70.1 ± 1.3 



LSS data added, the likelihood is quite symmetric around 
the central value. 

With this caveat in mind, we derive the following 2a 
(95.4%) constraints: 

- 1.75 X IQ-^^yr-^ < G/G < 1.05 x lO^^V"^ (H) 
and 

- 0.083 < 5G/G < 0.095 (12) 

We also plot 2D contours of marginalized likelihood 
distributions of ( versus G/G and SG/G in Fig. [5] and 



Fig. [TO] respectively. As expected, the variation of gravi- 
tational constant is strongly correlated with the value of 
C in this model. 

Some previous constraints on these two variables to- 
gether with the result of the present paper are summa- 
rized in Table [Til We note that in order to obtain such a 
constraint, one has to make some assumptions, either in 
the underlying theoretical model, or in the way G varies. 
This is particularly true for the case of constraints de- 
rived from CMB and LSS, as the impact of varying G 
on these are multitude. For example, Ref. 70] modeled 
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FIG. 7: One dimensional marginalized likelihood distribu- 
tions of G/G. 




SG/G 

FIG. 8: One dimensional marginalized likelihood distribu- 
tions of 5G/G. 



the variations of G by some hypothetical functions, while 
the present paper assumed Brans-Dicke model. One has 
to be careful when comparing the different limits, as the 
assumptions made are often different. Nevertheless, from 
this table we can get a feeling of the current limits on the 
variation of gravitational constants. 

IV. CONCLUSION 

In this paper, we use the currently available 
CMB (WMAP five-year^, ACBAR 2007l6|L, CBI 
polarizationfei] and BOOMERanG 2003[6l|66n63|) and 
the LSS data (galaxy clustering power spectrum from 
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FIG. 9: 2D contours of the marginalized likelihood distribu- 
tion of C, against G/G. 



0.2- 




-0.01 -0.005 0.006 0.01 0.016 0.02 0.026 0.03 0.036 



FIG. 10: 2D contours of the marginalized likelihood distribu- 
tion of (, against 5G/G. 



SDSS DR4 LRG data[6i|) to constrain the Brans-Dicke 
theory. We use the covariant and gauge-invariant method 
developed in paper I to calculate the CMB angular power 
spectrum and LSS matter power spectrum. 

To explore the parameter space, we use the MCMC 
technique. We parameterize uj with a new parameter, 
C = ln(l/w -I- 1), in order to explore the likelihood dis- 
tribution of the Brans-Dicke parameter w in a continu- 
ous interval. This method of parametrization is approx- 
imately equivalent to C = when a; is a large num- 
ber. It allows consideration of negative lo value, and also 
there is no arbitrary upper limit on (due to numerical 
problem, one has to choose a lower limit for |a;|). We 
explore in the range C S [—0.014,0.039], corresponding 
to w e [-cx),-71] U [25, oo]. 

We found that while the CMB observation could con- 
strain models with positive for the present data set 
and best fit parameter values, there is some degeneracy 
at a; < 0. The LSS data could effectively remove this 
degeneracy. Finally, using the CMB and LSS data, we 
obtain 2a (95.5%) limit on ( as -0.00837 < C < 0.01018, 
corresponding to a; < —120.0 or lo > 97.8. These lim- 
its may appear as weaker than previous limit obtained 
by Ref. [56| , even though we used later data than them. 
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TABLE II: Constraints on the rate of variations of gravitational constant. The errors are la unless otherwise noted. 



Parameter 


Value 


Method 


Reference 


G/G 


2±7 


lunar laser ranging 


MuUer & Biskupek 2007[71] 


[10-"yr-i] 


0±4 


big bang nucleosynthesis 


Copi et al. 2004 [72] 
Bambi et al. 2005 [73,1 




n i t 

U ± io 


helioseismology 


^1.^+1-. ^+ t^i 1 fine [7/11 
Ouentner et al. iyyo 




-6 ±20 


neutron star mass 


Thorsett 1996 [75] 




20 ±40 


Viking lander ranging 


Hellings et al. 1983 [76] 




40 ±50 


binary pulsar 


Kaspi et al. 1994 [77] 




-96 ~ 81 (2cr) 


CMB (WMAP3) 


Chang & Chu 2007 [70] 




-17.5 ~ 10.5 (2a) 


CMB±LSS 


Wu & Chen 2009 (This paper) 



However, this difference is largely due to the different as- 
sumption made in the constraint. Particularly, we con- 
sider case of a; < which was not considered in Ref . [s^ . 
As expected, the current limit on ui derived from CMB 
and LSS data is much weaker than those derived from 
solar system tests. However, the large temporal and spa- 
tial range probed by these observations make it a useful 
complementary to the latter. 

To examine whether the gravitational coupling is really 
a constant we introduced two new derived parameters in 
our MCMC code, one is G/G, the rate of change of the 
gravitational "constant" G at present, and the other is 
6G/G, the integrated change of G since the epoch of re- 
combination. We obtain the 2a limit for these two vari- 
able as -1.75x IQ-iV"^ <G/G< 1.05xlO-iV"^ and 
-0.083 < SG/G < 0.095 respectively. These limits are 
still somewhat weaker than the solar systems, but again 
they probed larger scales. Especially for this test, the 
assumptions made in each technique could be quite dif- 
ferent, which one must bear in mind when making com- 
parisons. 

The Planck satellite |80(| is expect to begin operation 
and bring back even better CMB data. The SDSS-3 
BOSS survey [Sll, WiggleZ [78|, and the LAMOST sur- 
veys are expected to measure galaxy power spectrum 



at higher redshift and with better precision. We look for- 
ward to obtain more stringent constraints on the Brans- 
Dicke theory and other scalar-tensor gravity models in 
the near future. 
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